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1. INTRODUCTION 

The prediction of boundary layer separation continues to be an important 
basic research problem in fluid dynamics with many applications in aerodyna 
mics, propulsion, fluidics devices, etc. The numerous theoretical and experi- 
mental studies that have been carried out to date have established three 
Important features of this problem that cannot be accounted for in the classi- 
cal boundary layer theory approach. (1) Viscous— inviscid interaction: the 

streamwise pressure gradient is not prescribed but is an unknown determined 
by the boundary layer displacement thickness distribution solution, (2) .Up- 
stream influence: the local behavior at some station x is affected by 

conditions do\^nstream as well as those upstream, the mathematical character 
of the flow being of an elliptic nature rather than parabolic (this is true 
even for a supersonic inviscid flow) . (3) Regions of reversed flow and possi- 
ble reattachraent : these introduce the well-known possibility of a singularity 

where the wall shear vanishes plus the numerical difficulties associated with 
imbedded regions of reversed flow near the surface. To complement the various 
numerical methods that are currently being developed to treat these features 
(Ref. 1,2 ■ ) the present paper describes an approximate analytical approach to 

the boundary layer -separation problem for subsonic" laminar separation bubble-type 
flows , ^ 

Our analysis is based on a triple-deck flow model which is an extension 

5 6 

of the earlier methods of Von Karman and Millikan, Stratford, and Curie to 
include the aforementioned effects of vlscous-inviscid interaction, upstream 
influence, and flow reversal and reattachment. The boundary layer is split up 
into two apx>ropriately-matched regions: a thin viscous sublayer region of 

negligable inertia near the wall underneath a thick outer layer of nearly 



1. INTRODUCTION 


The prediction of boundary layer deceleration and separation in regions of 
adverse pressure gradient continues to be an unsolved problem of fundamental 
practical importance in aerodynamics. Separation, for example, not only signi- 
ficantly influences the local flow field but also can influence the overall 
forces and moments acting on a wing or body. Although significant progress has 
been made in the theory of non-separated laminar and turbulent flows, existing 
methods have a number of deficiences which introduce serious drawbacks, when 
separation occurs? (a)' the flow model breaks down under conditions of vanish— 
ingly“small v/all shear and reversed flow: (b) they neglect the important 

effect of viscous-inviscid interaction induced by the boundary layer displace- 
ment thj.ckness growth, including the attendant upstream influence; (c) very 
little has been done to treat the turbulent case; (d) the numerical approaches 
involved are usually very lengthy and expensive to run and totally impractical 
for use in advanced engineering design calculations and parametric studies. 

As a first step toward remedying these deficiences, a new approximate three- 
layered theoretical flow model of boundary laj^er separation including viscous- 
ii).v''iscid interaction was conceived for the case of subsonic t\i?o dimensional 
steady laminar flox^ (see Fig. 1). In this approach, the boundary layer is 
split up into two appropriately-matched regions: a thin viscous sublayer region 

(having negligable inertia) near the wall overlaid by a thick outer layer of 
nearly inviscid (but highly rotational) flow. Inviscid interaction is accounted 
for by coup3.ing the perturbed inviscid flow to the total displacement thickness 
growth of the boundary layer 6*(x) using a linearized source distribution rep- 
resentation of 6*(x). The purpose of the research project was to study the 
theoretical development and numerical implementation of this "triple deck" 
flow model concept for laminar boundary layer separation problems, 



with the ultimate follow-on goalvof extending it to the case of ful3.y developed 
two-dimensional turbulent flow separation. 

2. PRESENT STATUS OF WORK 

The proposed triple deck interactive model for laminar flow has been 
completely worked out and some promising preliminary numerical results obtained, 
including development of a small perturbation analysis verifying the concept of 
the suggested iterative calculation approach. These results were documented 
and presented in an AIAA paper last summer; a copy of this, which contains all 
the details of the model analysis^ is given in the Appendix . 

During the remainder of last summer, the numerical feasibility of the model 
for treating separated flow conditions was further established, albeit in a very 
inefficient form unsuited for routine applications. Subsequently, working at 
a very low level of effort c^v/ing to a lack of student programmer time, we have 
been slowly developing a far more efficient algorithm to implement the inter- 


active model. 
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SUBSONIC laminar BOUNDARY LAYER SEPARATION 

AND REATTACHMENT WITH VISCOUS-INVISCID INTERACTION* 

Abstract 

A theory is developed for aonslmilar laminar boundary layers with separa- 
tion, reversed flow and reattachraent including global subsonic viscous-inviscid 
interaction and upstream influence with application to separation bubble pro- 
blems. The approach is based on an approximate triple-deck flow analysis which 
provides a unified, non-singular analytical model over a wide range of both 
attached and reversed flow states. The displacement effect on the invlscid 
flow is represented by a thin airfoil integral with a leading edge correction 
by Lighthill’s rule. A closed form solution of this model is given which shows 
the essential features of the problem and supports the numerical approach. A 
global iteration calculation method is devised which is free from the usual 
Crocco-Lees critical points occuring in interaction problems. Applxcatxpn to 
the case of a linearly-decelerating basic inviscid flow is discussed including 
some preliminary numerical results. 
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Parameter related to Inner layer thickness 
Characteristic length 

2(p-p^)/pu^ = 

^ o o 

2x/pu 

e 

o 

Parameter in inner layer solution 

Static pressure 

jsu^ x/u, pu^ jC/p 
o o 

x/x 


B 

Velocity component along x 

Velocity component along y (normal to surface) 

Coordinates along and normal to surface, respectively with 
X = x/yC being non-dimensional. 


p, basic c 


Boundary layer thickness 
Local Blasius thickness = 5.2i^^(l-C 
Displacement thickness . 

1.73>^/(1-C 
Momentum thickness 
Coefficient of viscosity 
Stream function, location of streamline 
Density 

Kinematic viscosity 
Wall shear stress = v(3u/3y) 


a/4 ^ 1/2 

, . ) Re 

p, basic c 


w 


Local flat plate value = ,332pu ^ (1-C , . 

p, basic X 

o 


Subscripts 

basic prescribed non-interacted pressure distribution values 

B Local Blasius (flat plate) value 



Cut-off point of basic pressure gradient 

End of Interaction Region 

Location of inner-outer layer interface 

Effective flat plate origin for favorable initial pressure gradient 

Effective flat plate origin for. downstream of interaction 

Maximum u station 
■ e 

Start of adverse pressure gradient 
Separation point 
Start of interaction 


Conditions at wall 



inviscld (but highly rotational) flow. Inviscid interaction is accounted for 
by coupling the perturbed inviscid potential flow outside the boundary layer 
to the displacement thickness growth 6*(x). There are several reasons for our 
selection of this approach. First, in the absence of viscous interaction 
effects the two-deck boundary layer model has been found to give very accurate 
predictions of both laminar and turbulent boundary layer behavior in adverse 
pressure gradients up to separation for a wide variety of flows including those 
with surface mass transfer,^ Hence its extension to include interaction, up- 
stream influence and separated flow would provide a valuable and relatively 
Inexpensive tool for making engineering analyses of separating flow. Second, 

our three-layer model approach is strongly suggested by the recent investiga- 

7 

tion of Stewartson and Williams, , who show that the asymptotic structure of 
incipiently-separati.ng laminar boundary layer flow at high Reynolds numbers 
does indeed have such a triple-deck character. Third, the relative simplicity 
of the analysis provides a clear insight to all the essential features and 
difficulties of the separation problem and lends itself to very straight- , 
foreward and Inexpensive numerical implementation. 

In Section 2 the basic simplifying assumptions, the analytical formula- 
tion and the salient features of the model are described. It is shown that 
the model has a hereto! ore-unrecognized double-branched solution character 
that enables it to pass through separation and reattacliment without singular!-- 
ties and provide a unified description of a wide range of both attached and 
separated flow states. Section 3 outlines a closed form method of solution 
which displays all the essential features of the problem and which also sup- 
ports our subsequent numerical approach. In Section 4, a global iteration 
technique of applying this model to the numerical calculation of subsonic 
viscous-inviscid interaction problems in adverse pressure gradients including 
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separation is then discussed, including some preliminary results for separa- 
tion bubbles occuring in a linearly-decelerating basic inviscld flow» We 
conclude in Section 5 with a brief discussion of further improvements and 
extensions of the present theory. 


2. FORMULATION OF ANALYSIS 


Assumptions 

VJe first introduce a number of basic simplifying assumptions, a) As a 
convenient idealization to illustrate the essential features of our approach, 
the flow is assumed two-dimensional and steady (neither of which is strictly 
true in practice) ’ and incompressible, although the present theory can be 
readily extended to compressible subsonic transonic or supersonic flow in- 
cluding heat transfer, b) The usual high Reynolds number boundary layer-type 
approximations are adopted Including the neglect of second order curvature 
effects. While this is not strictly valid near separation and reattachment 
points, it does not give significant errors unless the "fine structure" of 
the flow around these points is of interest. Indeed, the recent studies of 
Lees and Kleinburg,^ Stewartson,^ and Werle and Davis^^ have all shown that 
the use of boundary layer equations to predict separating flows is acceptable 
provided the pressure distribution is determined by vlscous-lnviscid interac- 
tion and the separation bubble region is thin, c) The flow is assumed to be 
laminar which is an idealization to illustrate the basic features of the 
analytical method and of course is of practical interest for low Reynolds 
number cases. The resulting analysis appears readily extendable to the tur- 
bulent case with the use of a suitable eddy viscosity model, since the essential 
physical features of the approach are independent of the details of the 
viscosity law. d) The global viscous-inviscid interaction is described by 
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coupling the perturbed inviscid flow to the boundary layer by means of a 
small disturbance source distribution (*‘thin airfoil”) model of the displace- 
ment thickness effect, thereby including consideration of the upstream in- 
fluence of downstream events (such as reattachment)*, (e) The types of flow 
considered here are those in which the basic non-interacting pressure distribu- 
tion basic zero in the region^ 0< x < followed by * an arbitrary* h: 7 

.i 

but finite length of continuous pressure rise over x < x < x , see Fig. 1 

C "■ ” o 

(this zero pressure gradient condition for x < x^ will be subsequently relaxed 

as discussed below. In the constant C , , region downstream of this cut- 

p, basic 

off point x^, the boundary layer is assumed to approach a flat plate behavior 
appropriate to this new pressure level. Although in the present study we em- 
phasize the case of reattaching laminar separation bubbles, the general approach 
is also applicable to other types of downstream boundary conditions such as 
those pertaining to a wake. 

Flow Model-Governing Equations 

Following the earlier ideas of Von Karman and Millikan as improved by 
Stratford and Curie, and more recently justified by Stewartsons asymptotic 
analyses of viscous interaction effects, the flow into an adverse pressure 
gradient region can be divided into three essential layers (see Fig, 1). The 
viscous boundary layer region for x > is divided into two strata; a thin 
"local equilibrium" sublayer 0 < y < near the surface where, inertia effects 
are negligible, overlaid by a thick outer layer of essentially inviscid but 
rotational flow with a small viscous total head loss that is conveniently 


* A vortex sheet model could also be used; in fact the present boundary layer 
model can be mated to any suitable numerical potential flow program (such 
as Hess and Smith^2) ^ 
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approximated by a local flat plate value. The third layer is the disturbed 
outer inviscld flow which globally interacts with the boundary layer. Since 


the analytical formulation of this two-layered boundary layer model has already 


been discussed in detail by several authors,^’^ only a brief outline heed be given 
here with emphasis on the new features not considered heretofore. 

Regarding the thick nearly-inviscid rotational outer flow layer 7 ^ Yj 
the velocity along any streamline x, ip is approximated by the variable total 
head Bernouli equation 




2(P-P^) 


( 1 ) 


where is the local flat plate (Blasius) solution pertaining to 

B r 5 

which approximately accounts for the viscous total head loss along the stream- 
line. Successive differentiations of Eq. (1) with respect to tp and neglecting 


the lateral pressure gradient yields the two equations 



( 2 ) 



The flow in the thin inner layer 0 < y < y^ (x) has negligible inertia and 
lateral pressure gradient (p^ - p^) and following Curie may be described 
by the following velocity profile which satisfies the momentum Eqn. as 
y 0: 


u - T y/y + (dp/dx) y /2y f a(x) 


N, 

y /u 


(4) 
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where x = v(3u/8y) Is the wall shear stress, the function a (x) is related to 

w 

the unknown value of as determined by matching with the outer solution, 

and N Is an arbitrary parameter with N > 3. The two layers are matched by 

requiring that ij>, u, 3u/3y and 3 u/9y each be^' continuous at y y^ , these 

conditions being sufficient to determine x, a(x) , y^ and with N appearing 

as an arbitrary profile parameter (the value N = 3.043 recommended by Curle^_ ' 

for xinseparated non-interacting flows In adverse pressure gradients is used 

■- 5 . • ■ 

here without change). As in Stratford's original work, the inner layer is 
presumed to lie within the linear region of the Blasius profile; consequently 
the aforementioned inner-outer matching yields the following approximate 
analytical relationship between the skin-friction and surface pressure dis- 
tribution: 


x^C 


(dCp/dx)" = . 


0122 (1-T) 


(N +3) 


Q;}-!) ^ (K^+3) 
3N^ (lH-1) (N-2) ^ 


(5a) 


= .0104 (1-T)^(1+2.02T) 


for N = 3.04 


(5b) 


where T = with = .332 basic emphasized 

that this relation is perfectly general as regards the nature of the pressure 
field; it refers either to a prescribed or freely-interacting ^^(x). 

Although Curie has shown that the assumption == y/p in the inner 
layer limits the pressure coefficient to values less than about 0,11, both his 
results and our own calculations indicate that accurate predictions of separa- 
tion from Eq. (5) may still be obtained for values of pressure coefficient at 


separation as high as 0.25. Moreover, since viscous-inviscid interaction 
further reduces the pressure level near separation (see below) , this limitation 



on Eq. 5 is not an important one in the present study. A more slgniftcant 
limitation stems from the neglect of the inertia effect (flo« history) of the 
inner layer; Eq. (5) consequently preducts that t = wherever the local 
pressure gradient vanishes, which certainly yields inaccurate results for 
prescribed pressure distributions with maxima or minima or for viscous-inviscid 
interactions involving local overshoots in the pressure. 

The unknown inner layer thickness is given by the aforementioned matching 

as 



with a limiting value of zero as the pressure gradient vanishes. The corre- 
sponding total boundary layer displacement thickness is composed of contribu- 
tions from the inner and outer layer regions as follows: 

-i:(- 

where 6 is the effective edge of the boundary layer (Appendix A). The dominant 
contribution to the displacement thickness is from the inner layer. Using 
Eqs. (1) ,, (4),and (6), Eq. (7a) becomes . — 

As shown in Appendix A, this can be approximated by the following expression 
which proves useful later: 
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(7c) 


«* =; 6„ + s - 


B 


which correctly passes over to the result 6* 6^* in the absence of pressure 

gradient (y^ -> 0^ where 6^^ and 6^ are the local flat plate displacement thick- 


nesses 


The effect of the global interaction between the boundary layer and the 
external inviscid flow is treated by correcting the original potential flow 
for the total boundary layer displacement effect following the approach pre- 
viously used successfully by a number of investigators, in which the velocity 
field perturbation due to the displacement thickness is represented by an 
equivalent source-sink distribution. Thus, for example, using the small dis- 
turbance approximation to Bernouli*s equation and evaluating the source-sink 
distribution effect at y - 6* as recommended by Preston (who found that 
this .gave good agreement with experiment for airfoil calculations) , the in- 
teracting pressure field can be expressed as tasic ^ where 


AC = 
P 


-2 


/i^^c 7 T" 

TT p, basic 


(x-5)^ (OdE 

(x-U^+6 (x)^ 


(8a) 


is the local correction to the prescribed basic (non-interacting) pressure due 
to the entire displacement thickness with ? the dummy variable of integration. 
Note that the integral is not singular at the source point C = x , which is 
advantageous in the numerical work. As shown in Appendix B, examination of the 
limits and integration by parts enables Eq. 8a to be expressed for computational 
purposes In the form 
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where x 0 and are the effective start and end, respectively, of the 
i f 

viscous— laviscid interaction region with x^ >> x^* Alternatively, using the 
classical "thin airfoil integral" formulation of AC^ in which the source-sink 
effect is evaluated along y = 0 instead of y = 5* with tne resulting leading 
edge singularity eliminated by the so-called Lighthill rule correction factor , 
Eq. (8a) reduces to the following well-known singular integral equation 
requiring the use of the Cauchy principal value at x = 



2 

ir 



(d6*/6x)dg 

x-^ 


(8c) 


where r = = 1..5 Re^ ^ = 1.5e^. Excluding the small higher order 

region of size r^ near the leading edge, this integral has the important pro-, 

perty of vanishing identically for the flat plate parabolic distribution 6* — 

5 * /*j • hence from Eq* 7C we see that only the pressure gradient effect 

B 

on 6* effectively contributes to the viscous interaction-induced pressure field. 


I mportant Features Pertaining to Separation 

For illustrative purposes, we consider the case of linearly-retarded flow 

u /u =(l - x)[C . = (2-x)x] for which an exact solution of the boundary 

15 

layer equations has been obtained by Howarth. The resulting solution of 
Eq. (5) for the shear stress distribution along the wall is shown in Fig. 2. 

From this Figure we can perceive several important features of the present approach 


tThis result may be understood from the fact that there is no pressure gradient 
immediately downstream of the nose of a parabola in incompressible flow. 
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First, we see that’ the two-layer boundary layer model is In excellent 
agreement with the exact solution, giving an accurate account of the decrease 

t 

in skin friction up to separation Including the location of separation itself 

around x w.l2. Actually, the two-layer result continues slightly downstream 
sep 

of this X corresponding to a small reversed flow (further comment on this is 
sep • 

given below) . 

Second, it is observed in Fig. 2 that the solution is double-valued : for 

each value of x and C one obtains two values of shear stress, one positive 

V — 

and one negative, corresponding to attached and reversed flows, respectively 
(note that Eq. 5b is a quartic In T with either a pair of real and a pair of 
imaginary roots or two pairs of imaginary roots) , Moreover the upper and lower 
"branches" of this solution are seen to be continuously connected. This 
heretofore-overlooked double-valued nature of the two -layer boundary layer 
model is analogous to the well-known double-branched similarity solutions of the 
boundary layer equations for negative values of the pressure gradient para- 
meter. Now, such a feature endows the present analytical approach with the 
inherent capability of giving a unified continuously-connected description of 
both attached and full-separated flow states provided we "unhook" the pressure 
field in the skin friction-pressure relation from a prescribed explicit de- 
pendence on X and instead allow it to be coupled to the boundary layer by 
viscous-inviscid interaction.* Thus, by interpreting Eq. 5 or Fig. 2 as a t 
vs. locus (Instead of t vs. x) the physical solution can move along AB 
with an increasing pressure toward separation, then from the upper to the 
lovjer branch through separation into a separated flow region along BC with 

* In the same manner as Lees and Reeves employed self-similar boundary layer 
solutions in supersonic viscous-inviscid interaction problems by "unhooking" 
them from the pressure gradient parameter. 
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decreasing interaction pressure, then back along CB through rising pressure, 
toward reattachment and through reattachment followed by downstream post- 
reattachment relaxation with decreasing pressure along BA. 

Third, it is important to note that unlike previous solutions of the full 
boundary layer equations, the approximate two-layer boundary layer model does 
not give a singularity at the separation point; this Is evident in Fxg. 2, 
where the slope of the wall shear stress at separation is seen to 

be finite (this can be also be verified by direct differentiation of Eq. 5 
with respect to x).^ Interestingly enough, this result is in qualitative agree- 
ment with experimental data and the predictions of exact numerical solutions of 

18—20 

the Navier-Stokes equations in adverse pressure gradients , both of which 

21 

indicate that dr/dx is finite (though large) as t -> 0. Lee has shown in 
fact that the usual boundary layer solutions giving a separation point singu-. 
larity are not valid in the immediate vicinity of when this is corrected, 

a regular behavior through separation is obtained. Indeed, he further shows 
by an asymptotic analysis of the Navier-Stokes equations near that to the 

leading order of approximation the viscous flow near the surface is described 
by precisely the inner layer flow model used herein. Additional support for 
at least the approximate correctness of the present model on this point can be 
obtained by examining the zero streamline ip = 0 whose locus indicates where 
the flow separates from the surface: integrating the inner velocity profile 

Eq. (4) with respect to y and using the matching relations, this locus is found 
to be governed by = 0 for t > 0 and in a reversed flow region by the 

relation ' 


^ The behavior near separation does not depend sensitively on N; for example, 
the cube root of the right side of Eq . 5 evaluated at x - 0 (which is 
approximately proportional to for linearly-retarded flow) decreases 

by only 16% when N is changed froS 3 to 4. 
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dC 


y “ 


6y 


N-1 


2 N-2 

m 


= 0 


(9) 


J l |;-0 


Differentiating this with respect to x and setting - 0 at t = 0, Eq. (9) 

predicts that the separation streamline leaves the surface at a finite angle y 

given by 


tany = 


dx 


= ~3 


^sep 


dx / dx 
dp/dx 


( 10 ) 


sep 


which under the assumption dp/dy = 0 is in exact agreement with the, well-known 

22 

value derived from the full Navier-Stokes equations by Oswatttsch, Both separa- 
tion and reattachment occur in a non-singular manner aat a well-defined acute - 
angle to the surface. 

Aside from Fig. 2., a fourth noteworthy aspect of the present analysis is 
the relative structure of the flow in a separation region predicted by the inner 
layer flox-J model. Thus, in addition to the aforementioned zero streamline, the 
zero velocity locus that delineates the extent of the reversed flow region is : 
found Eq. (4) to be governed by 


2Tf, 


f,B 


dC 

dx 


N-1 


y - 


N(N-l) N-2 


= 0 


( 11 ) 


u=0 


In addition the conesponding zero velocity gradient locus (3u/3y = 0) where the 
maximum reversed flow speed occurs is given by 


Tp 


f,B 


dC 

4 . L 

+ dx 


y - (N-l)« 


N-1 1 
r 

N-2 


= 0 


3u 


= 0 


( 12 ) 
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Since these locil derive from the inner layer solution, they must all lie 
within 0 < y < y^ and this is indeed readily verified to be the case by com- 
paring Eq. (6), Moreover, Eqs, (19), (11) and (12) predict that the ordinates 
(y)^ 0> (y) _Q and approximately in the ratio 3:2:1 independent 

of i^j this agrees with the relative size of these ordinates shown in Leal’s 
, 20 

Navier-Stokes solutions* 

Fifth, it should be noted that governing model equations automatically 
yield the correct downstream solution [namely a flat. plate boundary layer 
corresponding to C , . (x )] far downstream where the pressure gradient 

p ^ D€lSXC- 

vanishes; for when dC /dx 0 with x » x , Eqs. 5 and 6 give T ^ 1, y 0 

p c j 

while Eqs. (7) yields the correct corresponding displacement thickness. 

Sixth, shown in Appendix B the displacement thickness integral in Eq. 8b 
has the important property of becoming vanishing small (with a conccmlttantly— 
vanishing x-derivative) at both x = x^ and x = x^ provided and are 
chosen sufficiently small and large, respectively. Consequently, the inter- 
action-induced pressure AC^ and its streamwise gradient and hence (by Eqs. 5 
and 7) the corresponding shear and displacement perturbations all automatically 

vanish at x, and to any desired degree of accuracy, 
i It 

Allowance for an Initial Pressure Gradient 

The foregoing analysis can be extended to include boundary layer flows 

with a non-zero upstream favorable pressure gradient history in the region 

0 < X < x^ by adapting the equivalent flat plate initial condition technique 

of Curie. ^ Thus if the inviscid flow velocity distribution n^(x) increases 

to a maximum u at the pressure minimum station x = x , the boundary layer 
e m 

m 

velocity profile at x^ will be very similar to a Blasius profile based on 
some equivalent length ^ equivalent orig5.n is determined by 
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2 — ^ “ 

requiring the same momentum thickness 0* at x » Equating 0* ~ J 

^ o ui ' t!l O 

23 

u given by the accurate Tliwaites method to the equivalent flat plate 

value .45 (x -x ) yields 

m m L 



The subsquent boundary layer development in the adverse pressure gradient 
region may then be treated by the foregoing model provided x is replaced by 
X = X - x^ and the skin friction and the pressure coefficient are defined in 

terms of conditions at x^ instead of x^. 

The above equivalent origin approach can also be used to partially account 
for the "memory” effect of the upstream adverse pressure gradient history on 
the far downstream post-reattachment conditions in the case of the viscous- 
inviscid interaction separation bubble problem, as discussed below. 

3. ANALYTICAL SOLUTION 

An analytical solution of the foregoing theoretical flow model can be 
obtained which is a very helpful guide in constructing our numerical approach. 
To bring out the essential ideas involved, we adopt the approximate displace- 
ment thickness expression (7c) and neglect the 6* term in the integrand 
demorainator of Eq. (8) [neither of these simplifications alters the basic 
correctness of the subsequent analysis] ; then using the values of d and 

1/2 

given by Appendix A and Eq. (6), respectively, and noting that 6^* - 1.73 x e 

/. r n- 1/4 we obtain from Eqs. (5b), (7c) and (8) the 

' p, basic 

following trio of equations that characterize our interacting triple-deck flow 
model : 
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( 14 ) 




p ^ P „ “ /v-c 

p p , basic IT p, basic J 


(d6*/dx)dg 


(15) 


6* c: e 




(1-C ) 
P 


IM 


' «-4p\ .086(1-T)^(1-C 

•{1.73 +1 -rl) + — 

\ J (x dc Unr 


> = £ d(x) (16) 


Eq. (14) characterises the overall pressure-shear force balance of the boundary 
layer while Eq. (15) introduces the global effect of the displacement thickness 
growth on the pressure distribution in subsonic flow, as well as containing the 
"basic" non-interacting pressure distribution that characterizes the particular 

problem at hand. 

The trio (14) -(16) in general constitutes a difficult integro-dif ferential 
equation system to be solved for the three unkno\<mG T, and 6*. However, the 
presence of the parameter e and the fact that it is very small under precisely 


the high Reynolds number conditions appropriate to the present model suggests 
that a perturbation approach is feasible^ in which the solution is expanded in 
ascending powers of e, as follows: 


C = C , 4 + ^ 6C , + € AC « +. . . . 

p p, basic pi 


(17a) 


T = ■'basic + ^ ••• 


(17b) 


6* 


£ d, . + e d- + . . . 

basic 1 


(17c) 
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t A similar approach was successf uly used by Ting to solve a massive blowing 
problem having an analogous mathematical structure to the present one. 
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That :is , since the displacement effect' (within the accuracy of the boundary 
layer approximation) and hence the induced pressure are everywhere of small 
order e, the leading approximation is the non~inter acting solution; the small 


(order e) first interactive "correction” to this is obtained in terms of the 
6* for the "basic" flow. Thus, for example, substitution of expansions (17) 
into Eqs. (14) -(17) yield to first order that 


AC 


pl 




p, basic 


J o 


• X - ^ 




(18) 


(1-T^ . ) (1+2.02 T, . ) 

^ basic basic 

: - 3 - .2.02(1-4 T- . ) 
r'^basic 


AC 


£l 


“p, basic 


dAC _/dx 

+ 2 Ei 


dC I . /dx 
p, basic 


(19) 


with analogous expressions for AC^^ (iiivolving an integral of d d^/dx) and T^v 
etc.^ This procedure has thur converted the original problem involving an 


integral equation into a succession of quadratures of known displacement 
thickness functions. 

It is pointed out that the "basic" solution appearing in the foregoing 
analysis can be approximated by any reasonable distribution without 
altering the essential correctness or accuracy of the analysis, at least to 
first order in the interaction effects. Consequently, in the case where separa 
tion would occur in the basic flow, we may use an intelligent non-singular pro^ 
jective estimate of the 6* distribution across the length of the reversed flow 
in evaluating Eqs. (18) and (19). In this manner, our analytical solution 
provides the suggestion of (and support for) the streamwlse-pass iterative 
method of numerical solution discussed in the next section. 


— Xt is understood that e In e type terms may intervene between e and e ' In 
Eqs. (17). 


16 



4, NUMERICAL SOLUTION SCHEME 
Global Iteration Approach 

We have seen in Section 2 that the present triple-deck flow model has the 
inherent capability of giving a unified, continuously-connected description of 
both attached and fully-separated flow states provided we "unhook" the pressure 
field in the skin friction-pressure relation from a prescribed explicit der o 
pendence on x and Instead allow It to be coupled to the boundary layer by global 
viscous-inviscid interaction. The resulting solution can be thought of as 
being defined by the intersection of two locii (see Eig. 3) : the locus of all 

possible shear vs. pressure gradient values for the boundary layer as defined 
by Eq. 5, and the "interaction" locii of T vs. pressure at various x implied by 
the integral of 6*(T,Cp) in Eq. 8. The numerical solution of this problem is a 
difficult one since inclusion of the global subsonic interaction makes the pro- 
blem quasi-elliptlc and hence introduces the need to impose downstream boundary 
conditions. However, the discussion in Section 3 suggests that the most ex- 
pedient approach is to use a global Iteration method involving successively- 
refined streamwise passes. A very similar type of scheme has been proposed 
recently by Lees and Su for solving a stratified flow separation problem^^ and 
by Jobe and Burggraf for interacting trailing edge flow.^^ 

The manner of implementing this approach is as follows. Starting with a 
prescribed basic pressure distribution (which defines the specific problem), 
the complete non-inter active solution is calculated along the body, stopping 

ixi the event of incipient separation at some station x = x^ slightly upstream 

27 

of the separation point. Following Catherall and Mangier we then project 
from this station a first estimate of the dovTnstream behavior of 6* throughout 
the entire Interaction J5one assuming a regular continuous behavior, specifically. 
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we used a cubic polynorainal projection 6* = 6* ^ 

a (x“X- a = {d6^/dx) wtiile 3 .^ and a~ are determined by recjuiring 
3 1 1 ^ >5 

that 6* and d6*/dx be equal to local flat plate values at some first guess of 
the downstream interaction-termination distance x^. Having such a reasonably— 

J. 

behaved first approximation for 6*(x) over the entire flow length, one can 
make a corresponding smoothed first estimate of the displacement effect 
on the pressure distribution and thereby provide a basis for performing another 
streaniwlse pass in which the viscoas-inviscid interaction is now included. 
According to our analytical solution (Eq« 17) and previous discussion, the 
results so obtained for T. and ~ basic ^ fact be a fairly 

good first approximation to the final correct interactive values when the 
Reynolds number is large, provided x^ has been chosen sufficiently large. 

The second streamwise pass through is begun at some suitably small value 
of the upstream interaction-starting point x^. Retaining the aforementioned 
first guess for the numerically-smoothed , interaction-modified pressure dis- 
tribution C , = C ^ j + (AC ), (where subscript 1 denotes the first itera- 
pl p, basic p i 

tion) is used in Eq. (5b) to obtain an improved estimate of the skin friction 
distribution. Separation and penetration into the reversed flow region may now 
be allowed along the negative x branch of the x vs. curve (such as segment 
BC in Fig. 2); as the flow proceeds downstream the effect of interaction dies 
down and the solution moves hack through reattachment, relaxing toward a local 
flat plate behavior at the constant basic until finally the Interaction 
effectively ceases at the downstream location x^. Correspondingly, a nev7 cor- 
rected displacement thickness distribution can be calculated directly from 

t 

Eq. 7 (the polynomial projection now being discarded) and hence an improved 
interaction pressure correction (AC )^ for the second streamwise pass, and so 
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on. The modified pressure distribution for this second and all subsequent 
passes may be calculated by correcting the previous result, avoiding further 
direct use of C . . (x) ; thus for the n-th iteration with n > 2, = 

P y D3-S1C 

r /fc') + fAC -AC ) . These iterative streamwise passes are carried 

^ p,n p,n-l 

out until the solution no longer changes to within a specified amount. In 
general, the resulting skin friction and local pressure may not agree suf- 
ficiently with the desired downstream boundary conditions appropriate to the 
particular problem, implying that the original estimate of was too small. 

A series of larger values are then used, each accompanied by the aforemen- 
tioned sequence of streamwise pass calculations, until the downstream boundary 
conditions of vanishing AC^ and dC^/dx are satisfied to within the desxred 
degree of accuracy. 

There are several aspects of the aforementioned numerical scheme that 

should he noted. (a) Ck^ing to T and AC being defined in reference to the 

loca l basic flow, and the previously-noted property of Eq. 8 that both AC^ and 

dAC /dx can be made vanishingly small at a large and small enough x - and 
P 

X., respectively, it is seen that each streamwise pass inherently satisfies 
the correct initial and downstream conditions. In other words, our formulation 
insures that every iteration is automatically ’’tied down” to the proper non- 
interactive end points that define the basic problem, (b) In view of the 
experience of Jobe and Burgraf and others, convergence of the aforemen- 
tioned iteration process may be difficult to achieve x<rithout the use of ”under- 
relaxation,” in which only a fraction (usually half or less) o£ the correction 
to 6* is fed hack to obtain a new interacted pressure. A schematic illustra- 
tion of the proposed global iteration procedure involving 80% under-relaxation 
is shown in Fig. 4. (c) It is important to note that the present solution 


19 - • 



method does not encounter any Crocco-Lees throat-type singularity. This is 
due to the fact that we do not handle the viscous— inviscid interaction as an 
initial value problem (which Garvine has shown always leads to such a down- 
stream singularity) but instead as an effective boundary value problem: our 
n-th Iterative calculation of the interacting flow is always based on the 

known 6*(x) and C (x) distributions throughout the flow pertaining to the pre 
P 

vious n-1 streamwise pass, with the first streamwise projection done in the 


manner of Catherall and Mangier. 
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Shainroth and McDonald 
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have also 


shown that 


this approach provides an adequate treatment of upstream influence effects and 
downstream boundary conditions vjithout introducing streamwise saodle point-type 
singularities. (d) In carrying out ther solution, it is possible to further 
improve our treatment of the post-reattachment flow by incorporating a correction 


to the effective origin seen by the flow due to its* upstream adverse pressure 


gradient history. Thus, while in the case of an initial favorable pressure 


gradient the flow quickly forgets its history, the thickening and separation 
of the boundary layer associated with the traversal of an unfavorable^ pressure 


gradient region is not so readily forgotten even by the flow well dovmstream of 
reattachment. Hence the velocity profile at the interaction-cessation point x^ 
will have a Blasius shape but with flat plate displacement thickness ane skin 
friction values based not on but instead on x^ + where x^ is some effec- 
tive origin shift accounting for the upstream adverse pressure gradient history. 
Owing to the neglect of pressure gradient effects on the head loss in the outer 
layer (Eq. 1), the present theory does not include this memory effect. However, 
X can be estimated from the equivalent momentum thickness technique described 
in Section 2. Thus we obtain 
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u 

e 


5 


( 20 ) 






dx 


u 


e,o 


where as a first approximation in evaluating the integral, the viscous-inviscid 

interaction effect on the u (x) history has been neglected, i.e., we have used 

G- 

the given basic inviscid flow solution for u^» 


Some Frcliminary Results 

An attractive basic flow for theortical studies of separation is the 
special case of Howarth’s linearly-decelerating flow^"’ in which the basic adverse 
pressure gradient is applied right from the leading edge [that is, there is no 
boundary layer proceeding the pressure rise: = 0, 6*(x^) =0]. By cutting 

off this basic pressure gradient, reattachment and various finite lengths of 
separation bubble flow can be created (Briley^^ has obtained some exact numeri- 
cal solutions of the Navier-Stokes equations for several such cases which we 
also shall consider) . The downstream boundary conditions appropriate to this 
class of flows are that the flow attain an appropriate loc al zero pressure 
gradient with Blasius skin friction behavior at some distance downstream 
following reattachment where interaction ceases. 

The present investigation has been addressed to the specific examples of 

19 

linearly-decelerating inviscid flow considered by Briley invo3.ving tvjo dif- 
ferent lengths of basic adverse pressure gradient, one of which causes separa- 
tion and one of which does not. His skin friction predictions for these cases 
(labeled Two and One, respectively) are illustrated in Fig. 5. It is interesting 
to note that whereas Howarth's boundary layer solution predicts separation for 
both cases, Brileys results predict separation further dovmstream with in fact 
no separation at all in Case 1. This serves to emphasize that predictions of 
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separation by classical boundary layer solutions neglecting viscous-lnvlscrd 
interaction are not always reliable. Note also from Case 2 that both separa- 
tion and reattachment occur without any singularities. 

An illustration of the typical displacement thickness obtained by the 
polynomial-projection procedure used in the first streamwise pass of our pro 
posed iteration method is given in Fig. 6 for Case 2. Our first streamwise- 
projected 6* distributions for both Cases 1 and 2 (based on the crudefirst 
estimate x, - .50 - see Fig. 5) are compared with Briley's results xn Fxg. 7. 

In general, it appears that this "first guess" procedure yields reasonably 
good results with a reasonable choice of x^. Note that the present theory does 
not give quite the same value of the downstream local Blaslvs thickness because 
of its very approximate account of the upstream history effect on the effectxve 

origin. 

The first interactive pressure corirection AC^ associated with the integral 

of our project 6* for Case 2 is &ho\m. in Fig. 8a. This is a rough result only, 

having been obtained with only a crude numerical integration routine and devoid 

of the Lighthill correction for the leading edge singularity. Nevertheless, 

when a smoothed-out curve fit of this result (shown dashed) is subsequently 

added to C . , the resulting first interactive pressure distribution C - 

p, basic 

(Fig. 8b) is quite reasonable: the effect of interaction has appreciably 

smoothed and broadened out the basic pressure distribution so as to relieve 
the adverse pressure gradient preceeding separation (thereby explaining why 
separation tends to be significantly delayed compared to classical boundary 
layer theory predictions). It can also be seen that this interacted pressure 
distribution tends to form a plateau in Che separated flow region, although 
this tendency is not as pronounced here as it would be in supersonic flew owing 
to the difference in shapes that constant pressure surfaces assume in subsonic 
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vs, supersonic irrotational inviscid flow. Downstream, reattachment occurs- 

through a continuous taonotonic pressure rise, A13. these features are in 

32 

general agreement with experimental observation. 

Our current efforts are involved with improving the original implementive 
techniques of the global iteration procedure (as regards a more accurate numeri— 
cal evaluation of the interaction integral and incorporating under relaxation 
into the iteration scheme) so as to extend the foregoing preliminary results 
to iteratively-converged final answers including wall shear distributions and 
velocity profiles along the entire length of the interaction and separated flow 
regions. Upon achieving this, it is then planned to run some parametric studies 
of the linearly-decelerating problem varying the length and shape of the adverse 
pressure gradient region, the value of basic^^c^’ Reynolds number, and 
allovjing non-zero initial boundary layer thickness (x^ > 0). Hopefully, some 
cases might also he run where at least limited comparisons could be made with 

33 

experimental data on airfoil leading edge laminar separation bubbles 

5 , CONCLUSION 

Although approximate the present theory has the virtues of sound physical 
modeling of the essential flow features including viscous-inviscid interaction 
and upstream influence, an analytical formulation which is readily implemented 
numerically, and good engineering accuracy. The triple-deck model has the in- 
herent capability of passing smoothly through a separation point into a reversed 
flow region and back through a point of reattachment without singularities. In 
addition to providing a method for obtaining approximate engineering solutions 
of boundary layer separation problems, the present theory would bje a useful 
analytical tool for simulating the process of reattachment itself by providing 
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suitable fully-separated shear flov configurations. Furthermore, the present 
theory may serve as a valuable aid in the development of more rigorous and 
detailed exact numerical methods for high Reynolds number flows by providing 
estimates of the location and mesh size requirements of the high shear and 
reversed flow regions. Finally, we reemphasize that the triple-deck model 
approach is applicable not only to flows with reattachment but also appears 
adaptable to those where other types of downstream condition pertain such as a 
wake . 

Many improvements and extensions on the present analysis appear possible. 
For example, the assumption that the total head loss in the outer layer is not 
influenced by the pressure gradient can be relaxed, as can the neglect of the 
convective inertia effects in the inner layer (this would enable the theory to 
treat basic pressure distributions V7ith maxima or minima) . Extension can also 
be made to include' compressibility effects for either subsonic, transonic or 
supersonic flow, including heat and mass transfer. Furthermore, it xs clear 
from Stratford’s pioneering study of the two-layer model for non-interacting 
incipiently-separated turbulent boundary layers^ ^ that the present approach can 
extended to interacting turbulent flows provided a suitable eddy viscosity 
model for the strongly-adverse pressure gradient and reversed flow regions is 
available. Finally, the theory ban be applied to treat three dimensional flows 
since there is nothing inherent in the character of either the inner or outer 
layer approximations that cannot be extended to include the presence of a cross 

flow. 
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APPENDIX A 


Approximate Expression for Displacement Thlcloiess 

A simplified version of Eq. (7b) can be derived by developing an analyti- 
cal approximation to the integral term as follows. Noting that the radical 
in the integrand is simply that the equality of stream function 

requires u ^ dy = obtain 

* outer B B 



(A-2) 


where n = y/^ , 6 and 6^* pertain to the local basic inviscid flow u^(not 

u ), and the Karman-Pohlhausen polynomial approximation - 2ti - 2n + 

eo fcj 

has been used to evaluate the integral. Then by neglecting 3ud as 
small compared to unity, wtd obtain upon substitution of A-2 in Eq, (7b) the , 
final result given by Eq. (7c). 

In the aforementioned expressions for 6*, the value of the conespandlng 
boundary layer thickness 6 can be estimated from the well-known Thwaites 
formula as 



o 




(A-3) 



APPENDIX B 


Properties of the Displacement Thickness Integral 

The induced pressure field is governed by the global effect of the dis- 
placement thickness as shovm in Eq. 8* Assuming for the moment that the viscous 
interaction-induced pressure gradient effect on 6* is confined to some region 
^ X < the Integral in this Eq. can be written as 




0 (x-^) -HS*" 


J 


d6 * 




Xj[ 


<» . 
Xf 


(B-1) 


where x. (->0) and x » x are the effective start and end of the interaction 
i X c I 

region, respectively. Here, the first and third integrals account for the con- 
stant pressure regions of the basic flow while the second is the contribution 
of the interaction region. We nov^^ consider Eq. B-1 applied to some stream- 
wise point X within the interval 5 x < x^ and proceed to examine each of the 
three integrals on the RHS. 

Taking note of the fact that 6* ^ 1.73e (x v'l^C 7 and choosing 

B ' P, basic ^ 

^i enough that ^^i^ ^0, the first integral can be rewritten 

as 


•-.X, 

^ / 1 
'2 

1.73e 

J o 


(x-eUd6jj*/dx)d^ 

(x-e)^+6*^(x) 
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/?[(x-5)^+6*^(x)] 


(B-2) 


Now since <S*"e is small compared to unity and since Eq. B-1 pertains to 
X > x^, 6* makes a significant contribution to the integrand only within a , 
distanceJ^^5^6^*(x^) of the end y^oint ^ - x^; consequently , this RHS integral can 
be further decomposed into the two parts 
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1 
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dC 


(x-C) 


A7 


^*1 


j 


X ‘ 


(x-o dc 

2 2 


which standard integral tables evaluate as 



x+(x^-Jl)+2/x(x^-f) 

1^ 1 w 

f* ' ? 2 


“< log 

2^ 1 

(x-x.)^+il^+25l(x-x.) 

j, 1 
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We note that the second term vanishes while the first gives a non-singular 

result at X -> x^» Most importantly, we note that choosing = A == 

eliminates the first term at all x > x^, leaving only the second term which 

vanishes at x = and as x >> x^. • The first integral in B—1 can therefore 

t 

be made negligiable by this suitably small choice of x^ . 

Substituting the value of <5^^ into the third integral on the RHS of (B~l) 
and proceeding in the same manner as we did vrith (B-2) , this integral can be 
written 


2 

fe-o to 1 

oo 

r d? , 1 

(x“^)dC 

1,73e 

J 

(x-D^+6*^- J 

i 

x^+k (x-C) 

2 2 
x^. (x-O + ^ 
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where the 6*'^* term in the. integrand is presumed to have a negligable effect 
beyond some small distance k from x^. Again using standard integral tables, 
the RHS becomes 


log! 


x^+k-x 

f 


/ x^+k+x-2 /x ( x^+k) j 


f(x.+k-x)Vk^ 

+ log 2 ^2"^ 

2v^ I +k 


^This choice of x. is analogous to the use of a Lighthill correction factor in 
Eq, -8c to wipe out the leading edge contribution in the alternative singular 
airfoil integral formulation. 



which we see vanishes at x — and as x << x^ for all k provided x^ is taken 
sufficiently large. 

We have thus shown that provided and are chosen sufficiently small 

and large, respectively, only the second integral on the KHS of Eq. B>-1 does in 

fact effectively contribute to the induced pressure field as assumed, and that 

ACp vanishes at these end points. Moreever, since this type of integral is in 

effect a solution to LaPlaces Equation, dAC /dx inherently vanishes at these 
- '26 ^ 

pornts as v/ell. In carrying out numerical evaluations of the second integral, 
it proves convenient to elimraate the derivative of S* through integration by 
parts so as to obtain 


(x-C)^ d? 
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[(x-C)^-6*^16* 
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,;Fig. 4, Global Iteration. Procedure. 
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